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We report the study of a new experimental granular Brownian motor, inspired to the one published 
in [Phys. Rev. Lett. 104, 248001 (2010)], but different in some ingredients. As in that previous work, 
the motor is constituted by a rotating pawl whose surfaces break the rotation-inversion symmetry 
through alternated patches of different inelasticity, immersed in a gas of granular particles. The 
main novelty of our experimental setup is in the orientation of the main axis, which is parallel 
to the (vertical) direction of shaking of the granular fluid, guaranteeing an isotropic distribution 
for the velocities of colliding grains, characterized by a variance Vq. We also keep the granular 
system diluted, in order to compare with Boltzmann-equation-based kinetic theory. In agreement 
with theory, we observe for the first time the crucial role of Coulomb friction which induces two 
main regimes: (i) rare collisions (RC), with an average drift (a;) ~ Vq, and (ii) frequent collisions 
(FC), with (ijj) ~ vq. We also study the fluctuations of the angle spanned in a large time interval, 
A6, which in the FC regime is proportional to the work done upon the motor. We observe that 
the Fluctuation Relation is satisfied with a slope which weakly depends on the relative collision 
frequency. 



PACS numbers: 45.70.-n,05.40.-a 



I. INTRODUCTION 



Brownian motors (BM) are devices that can rectify 
thermal fluctuations, in order to perform work against 
external loads [l|, ^ . The basic underlying mechanism re- 
lies on the presence of non-equilibrium conditions, break- 
ing the time-reversal symmetry in the dynamics, together 
with some spatial anisotropy, which allows unidirectional 
motion. Although these general constraints are clearly 
understood [3|], many open questions remain to be an- 
swered, concerning the several different mechanisms for 
the realization of such devices. For instance, given a par- 
ticular shape of the probe, the prediction of the drift di- 
rection is far from obvious, in particular if several sources 
of dissipation arc present in the system, inducing com- 
petitive effects. 

Recent years have seen an increasing wide interest on 
Brownian motors directly inspired to the original setup of 
the Feynman's ratchet [3[. In these "coUisional Brownian 
motors" (CBM) fluctuations are induced by the collisions 
of an asymmetric probe with particles of molecular flu- 



ids at different temperatures [3-|g|- Since the presence 
of dissipation is a fundamental ingredient to induce non- 
equilibrium conditions, a natural framework where these 
kinds of systems have been studied is the realm of granu- 
lar media, where interactions do not conserve energy due 
to inelasticity^ Several experimental |7Hl0l| and theoreti- 
cal results pTl - ll4l | have been obtained for these systems. 

More recently, another source of dissipation has been 
shown to play an important role in the dynamics of CBM: 
the Coulomb (or dry) friction |15l - [l7| . Its main effect is 
the introduction of two dynamical regimes, where the be- 
havior of the systems is dominated by collisions or fric- 
tion, respectively. More surprisingly, it has also been 
shown that the Coulomb friction itself can be sufficient 
to drive a motor effect, even if the probe is in contact 
with a single molecular fluid at equilibrium [o, lla] • The 
role of friction has also been studied in other systems 
showing motor effects [l^ [20] , where fluctuations are in- 
troduced by noise terms, which are not related to particle 
collisions. 

The context of granular systems also paves the way to 
the realization of experiments aimed at validating some 
important general relations derived for nonequilibrium 
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systems, such as the Fluctuation Relation [21| or the 
Hatano-Sasa relation [22| . In particular, in granular sys- 
tems, where noise and time-scale separation are often not 
fully under control and where some coarse-grained on the 



accessible quantities is present, the experimental study of 
these relations is very useful to assess such results in more 
general situations |23l - [25| . 

In this paper we consider a new experimental setup for 
a frictional granular CBM, in order to get closer to condi- 
tions where kinetic theory can be applied. Moreover, at 
variance with previous studies, we also take into account 
the presence of Coulomb friction which induces interest- 
ing behaviors. This allows us to compare experimental 
results with analytical predictions of kinetic theory. Fur- 
thermore, we focus on the study of the non-equilibrium 
fluctuations of the spanned angle in a time interval, which 
is related to the work done by the CBM. Our findings 
suggest that a symmetry relation for these fluctuations 
is verified in our system, in agreement with previous re- 
sults for similar experiments |23j . 



II. SETUP 

The two main components of our setup are the gran- 
ular gas and the rotor, see Fig. [T] for visual explanation. 
The granular gas is made of A^ = 50 spheres of poly- 
oxymethylene (diameter d = 6 mm and mass ?7i = 0.15 g) 
contained in a polymethyl-methacrylate (PMMA) cylin- 
der of circular base, with area A « 6.36 x 10'^ mm^ and 
maximum height 70 mm. The cylinder is shaken with a 
sinusoidal signal at 53 Hz and variable amplitude, which 
is measured by the maximum acceleration F = a^ax/G 
rescaled by the gravity acceleration G. The velocity dis- 
tribution of the spheres on the plane perpendicular to 
the rotation axis is obtained by particle tracking via a 
fast camera (see [2y| for details on the procedure) and is 
fairly approximated by a Gaussian, 



Pai^) 



-K/i^<) 



(1) 



where z is the coordinate of the 



where the "thermal" velocity vq has been introduced and 
Vx may be replaced by Vy because the system is isotropic 
on the xy plane. Small deviations from the Gaussian 
are observed but are neglected for the purpose of this 
study; see [l3| for details. We have changed F from 5 
to 21, finding for vq values from 120 mm^ s~^ to 500 
mm^ s^^. The vertical average profile of the gas de nsity 
is close to the so-called "leading order" distribution [231 ■ 

vertical axis. This is sufficient to evaluate the average 
density surrounding the rotator, estimated to be n = 
77.0(2:*) where z* is the mid-heigth of the rotator. 

Suspended into the gas, we have a pawl (also called 
"rotator") rotating around a vertical axis. The rotator 
is a PMMA cylinder with a rectangular base (see lower- 
right inset of Fig. [T]) of perimeter 5", height h and total 
surface of the sides 'E ~ Sh ~ 1.2x 10^"^ mm^. The rota- 
tor (including its axis) has mass M = 6.49 g and moment 
of inertia / — 353 g mm^. The axis of the rotator is sus- 
pended to two spheres bearing. The position vs time 
of the rotator is recorded by an angular encoder (Avago 



Technologies, model AEDA-330), enclosing the bearings. 
In Fig. [TJ top-right inset, the definition of some quan- 
tities useful for the theory can be found. The study of 
the Brownian motor phenomenon is obtained by apply- 
ing insulating tape to the rotator, partially covering its 
two largest surfaces, see lower-right inset of Fig. [T] for an 
explanation. Turning upside-down the rotator allows one 
to invert its chirality, defining two possible orientations 
which we call "right" and "left" . 
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FIG. 1. Color online: the experimental setup 



A close analysis of the dynamics of the rotator shows 
that the angular velocity uj is well described by the fol- 
lowing equation of motion: 



w(t) = -Acr[a;(f)] - 7aw(t) -f- T]^oii(t) 



(2) 



where A = Ffrict/I = 38 ± 4 s^^ is the frictional force 
rescaled by inertia, 7^ = 6± 1 s^^ is some viscous damp- 
ing rate related perhaps to air or to other dissipations 
in the bearings, and r]coii{t) is the random force due to 
collisions with the granular gas particles. The granular 
gas itself is stationary and (roughly) homogeneous with 
a mean free path proportional to A = (nS)~^. The pawl 
is further characterized by its symmetric shape factor 
{g^)surf = 1.51 (see [131 for details), with Qsurf being a 
uniform average over the surface of the object parallel to 

the rotation axis and g{s) = ''^^j, ^^' with £(s) = £ x n{s) 
which is the unit vector tangent to the surface at the 
point r(s), and n{s) is the unit vector perpendicular to 
the surface at that point. We have also introduced the 
radius of inertia Rj = ^JljM of the rotator. We refer 
to the top-right inset of Fig. [T] for a visual explanation 
of symbols. The restitution coefficient between spheres 
and the pawl has been measured to be otj^ sa 0.67 on the 
PMMA naked surface and a_ w 0.35 on the tape-covered 
face. 

It is useful to introduce the "equipartition" angular ve- 
locity Wo = VQi/Ri where e = y/^- Note that, because 
of inelastic collisions and frictional dissipations, the rota- 
tor does not satisfy equipartition and wq is only a useful 



reference value. It is natural to adimensionalize the ro- 
tator angular velocity defining 51 = -^. 



been obtained, for a gene ric asymmetric rotator in a di- 
lute gran ular gas, in [l3|. For our particular shape they 



read 
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III. THEORY 

The single particle probability density function (pdf) 
p{uj,t) for the angular velocity of the rotator is fully de- 
scribed, under the assumption of diluteness which guar- 
antees Molecular Chaos, by the following linear Boltz- 
mann equation [IJ, [l5|, II4I 

dtpiuj, t) = d^ [( Acr(w) + '-fauj)p{uj, t)] + J[p, pg] (3a) 

J[p,Pg]^ fdLu'W{Lu\L0')p{L0\t)^piLU,t)f,iLu), (3b) 



(3c) 

(3d) 
(3e) 



fc{uj) = / duj'W{uj'\uj) 



ds 



Wiio'luj) = pS I — I dvp<,(v)e[(V(s) - v) • n]x 



|(V(s) - v) • h\S[uj' -uj- Aa;(s)], 
[Y{s) - v] • n g{sy 



Aa;(s) = (1 + a)- 



Ri 



l + e2g(s)2 



where wc introduce the rates Vl^(aj'|a;) for the transition 
UJ — )■ w', the velocity-dependent collision frequency fc(j-o), 
the pdf for the gas particle velocities P(,(v) and the so- 
called kinematic constraint in the form of Hcaviside step 
function 0[(V — v) • n] which enforces the kinematic con- 
dition necessary for impact. Here V(s) = uz x r(s) is 
the linear velocity of the rotator at the point of impact 
r(s). The collision rule is implemented by Eq. ([5e)) |13| . 
An estimate of the ratio between the stopping time due 
to dissipation (dominated by dry friction) t^ '^ ^ and 



the coUisional time Tc 



nEuo 



is given by the parameter 



n-i ^ e"Sw^ 



TA 



V2Tr Ri A Tc 



(4) 



This parameter controls the transition from a regime (at 
/3^^ <C 1) with fast stopping due to dissipation, called 
RC (rare collisions), and a regime (at f3~^ ^ 1) with 
the rotator always in motion, continuously perturbed by 
collisions, called FC (frequent collisions). 

When the mass of the rotator is large with respect to 
the mass of the granular gas particles, collisions are small 
perturbations to uj(t) (see Eq. (I3el)) . Then, it makes sense 
to expand Eq. (jH]) in powers of e [i2:ll3l: by retaining only 
up to the second derivative, a Fokker-Planck equation 
is obtained. The basic result of this procedure is that 
the coUisional noise in Eq. ([2]) is cast into the sum of 
a white noise T]{t) plus a viscous drag and a systematic 
force inducing the motor effect: 



Vcollit) -^ V{t) - lg^{t) + Tynoto 



(5) 



with (77) ^ and {r]{t)r]{t')) = TgS{t - t'). The expres- 
sion for 7g, T„iotor and the amplitude of the noise Tg have 



7s 



'^moto 



ISir 1 a_|_ — a_ 



32 e2 + a+ +a_' 



rg = (l + a)7g— z;2, 



(6) 
(7) 
(8) 



with a = (q;+ + a_)/2. We want to highlight, here, a 
fundamental point concerning the coUisional noise rjcoif- 
as Eq. ^ explicitly shows, such a noise is in general not 
white, and, even more importantly, it is not independent 
from the instantaneous velocity w. This makes sense, as 
it is the superposition of the variations of angular velocity 
due to collisions, which, as shown by Eq. (pc|) . depends 
on UJ. For this reason our model, described alternatively 
by Eq. ^ or Eq. ([3]), is very different from a model ~ 
apparently similar - recently introduced in [201, ^^ well 
as from other previous models |28j . 

In the FC limit, f3~^ ^ 1, the Coulomb friction term 
and the external viscosity may be neglected, i.e. 

7aW -I- Aa{uj) < jguj, (9) 

so that Eq. (21) is cast into the much simpler form 



CJ(t) = -"fgUj{t) + ri{t) + Tjnotor- 



(10) 



From such an equation one may estimate the average 
velocity of the Brownian motor to be jla, [Tg 



Tjnoto 
I9 



(11) 



Another interesting observation follows from Eq. ((TU)) . It 
concerns the fluctuations f{A6) of the spanned angle in 
a time interval of length At, AO = 9{t + At) - e{t) for 
any t in the steady state. For the particularly simple 
linear Langevin case, Eq. (|10p . it can be shown that such 
fluctuations obev for large At, the following Fluctuation 
Relation (FR) [2g: 



0(A0) = log 



/(Ag) 
/(-A0) 



sAO, 



with 






(12) 



(13) 



Wc mention that such an FR is closely related to the 
FR for the entropy produced in the time At, which in 
this system is approximated by the work done by the 
"Brownian motor force" W sa T,notorA9 divided by the 
"temperature" (w^) [3G |. 

In the RC regime, on the other side, one may assume 
that the probe's dynamics is a sequence of independent 



kicks received at zero velocity, resulting in the follow- 
ing formula for the adimensional average angular veloc- 
ity [H: 



-A r 
4 ^ 
/tan^i V? 



(14) 



g^^[(l + a+f-(l + a_)2] (15) 



V V? 1 + e 



(16) 



where £, = mL^/{4:I) and L is the lenght of the pawl. We 
note that for the dimensional angular velocity this means 
|a;| ~ VqC. In the RC regime the behavior of (/)(A6') is 
unknown in principle. A FR for the entropy production 
certainly exists, but we are not aware of a simple relation 
between A9 and the entropy produced in a given time 
interval. 



IV. MOTOR EFFECT 
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FIG. 2. Color online: trajectories with different shaking 
parameters and different chiralities. The inset shows the be- 
havior of the symmetric (i.e., without tape) pawl. 



In Figure [2] the evolution in time of the angle 9{t) 
spanned by the rotator is shown for different choices of 
the maximum shaking acceleration F and different ori- 
entation of the asymmetry ("left", L, or "right", R). A 
measure of the "thermal" velocity of the particles, vq, 
through the fast camera, allows one to determine (3~^ 
which estimates the relative relevance of collisions with 
respect to Coulomb friction. A steady drift, signaling 
the presence of Brownian motor effect, is observed both 
in the friction dominated regime (/3^^ < 1) and in the 
collisions dominated regime (/3~^ > 1). Turning the ro- 
tator upside down, i.e. changing its asymmetry from L 
to R or viceversa, inverts the sign of the drift: such an 
observation confirms that the observed drift is caused by 
the asymmetry of a, as expected from the given theoreti- 
cal arguments. A further confirmation that the observed 



drift is due to the surface heterogeneity comes from the 
study of the "symmetric" rotator, which has no patches 
of insulating tape: on the same timescale only a very 
weak drift is observed, much weaker than the one ob- 
served with the L- and R-type rotators. We impute such 
a weak bias measured with the symmetric rotator to im- 
perfections of the experimental setup: any kind of asym- 
metry (e.g. not sufficiently precise vertical alignments, 
not perfectly circular profile of the container, etc.) may 
induce secondary motor effects. 
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FIG. 3. Color online: Average angular velocity, rescaled by 
ojo, of the rotator: experiments, theory and simulations, for 
several values of /3~^ and both chiralities. In the inset only 
simulations are show, in order to appreciate the comparison 
with both RC and FC theoretical limits at very small and 
very large /?""'. 



In Figure [31 the rescaled average angular velocity of 
the rotator, (f2), measured in experiments with both L- 
and R-type rotators, is shown as a function of /3~^. The 
velocity of the L- rotator is changed of sign for the purpose 
of a better visualization. In the same plot we have also 
shown the results of numerical simulations of Eq. ([3]) with 
the same parameters as in the experiment (only for the 
R-type rotator). Details about the simulation can be 



found in the Supplemental Material of ref. [17|. As a 
first observation, we highlight the quite good agreement 
between experimental and numerical data. 

The thick blue dashed and magenta dot-dashed lines 
represent the theoretical predictions for the RC and FC 
limits, Eqs. ([T^ and PT|) respectively. In the inset of the 
Figure we have displayed the results of the simulations 
on a much wider range of /3~^, in order to appreciate 
the agreement with the theoretical limits. The simula- 
tions teach us that such theoretical predictions for the 
RC limit (FC limit) are useful for quite small (large) val- 
ues oi f3^^. The experimentally accessible values of /3^^ 
appear to be at the crossover between the two regimes: 
nevertheless they span a sufficiently wide range, so that 
both the (n) ~ /3~^ behavior (RC regime) and the trend 
toward saturation (uj) — >■ const. (FC regime) can be iden- 



tificd. We consider this to be the best comparison, up to 
our present knowledge, between experimental granular 
Brownian motors and kinetic theory. We mention that it 
is quite difficult to expand the accessible /3~^ range. In- 
deed exploring smaller values of /3~^ requires a consider- 
able increase of the dry friction coefficient A, which is not 
under our direct control; moreover, a large A may am- 
plify non-ideal effects where the behavior of the spheres 
bearings do not follow the Coulomb law: such effects arc 
already observed here (at small velocities small devia- 
tions are observed) and are likely responsible for the not 
perfect match with simulation results at small f3~^; we 
also notice that static friction is not considered here, but 
it could become relevant at large A [31|. The opposite 
limit, i.e. large values of I3~^, are even more difficult 
to be attained, since they would require a larger colli- 
sion frequency: however the maximum acceleration of 
the shaker is a hard limit, while increasing the number 
of grains does not trivially produces the desired result, 
for two reasons: 1) higher densities correspond to a re- 
duction of the average kinetic energy and 2) Molecular 
Chaos is only guaranteed at low density. Notwithstand- 
ing those limits, we believe that the results of Fig. |3]are 
already a quite good success of kinetic theory and make 
us claim that Eq. ([3]) is a fair description of the experi- 
mental setup. 



cases we have also superimposed a parabolic fit of the 
data at the largest available time. We mention that the 
characteristic time Tc for the decay of the angular ve- 
locity autocorrelation, not shown here, is in the range 
0.03 — 0.06 s. As frequently observed [S^, [3j1, the em- 
pirical large deviation rate tends to become independent 
from time only at very large At 3> Tc- Here we evaluate 
At ~ 1 s as a minimum value before considering reached 
the large deviation limit. In both experiments we may 
appreciate deviations from the parabolic fit, i.e. slightly 
non-Gaussian tails, signaling that we are indeed probing 
the large deviations of the pdf of A9. 
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V. NON-EQUILIBRIUM FLUCTUATIONS 



FIG. 5. Color online: Asymmetry function (l>{A9) for the 
pdfs of A9 for two experiments in the RC regime (left) and 
in the FC regime (right), at several values of At. 
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FIG. 4. Color online: Empirical large deviation function for 
the fluctuations of A9 for two experiments in the RC regime 
(left) and in the FC regime (right), at several values of At. A 
small correction is operated to each pdf f{A6), normalizing 
it by its maximum value, before taking the logarithm, for the 
purpose of a better vertical alignment. 



In Fig. |4] we have displayed the empirical large devi- 
ation rate (Idr) - log[/(A6')]/At of the pdf /{AO) for 
different choices of the time window At, in two experi- 
ments with a small and a large value of l3~^. In both 



The asymmetry function 0(A0) = log 



/(Afl) 



/(^AeJJ for the 
pdf of A9 is shown, for the same two experiments, in 
Fig. [5l At values of At large enough, but smaller than 
those necessary to achieve a stable large deviation rate, 
the asymmetry functions already display a linear behav- 
ior ^ sA9 with a slope s not dramatically changing with 
time. The values of the slope have been measured for all 
experiments and many values of At and are reported in 

Fig.m 

While confirming the direct observation done in Fig. [SJ 
Fig. |6] clearly teaches us that the slope of the (f){Ad) be- 
comes stable at times of the order At ~ 0.3 s. In the inset 
of Fig. [5] we have displayed the slope s (from a constant 
fit of the main plot including values At > 0.32 s) vs. 
P^-^. The squared symbols joined by the green dashed 
line represent the result of formula (|13p . Many impor- 
tant comments are in order here: 1) for large values of 
f3~^, as expected, dry friction becomes negligible and - at 
the leading order - the system is described by Ec^. (|10p . 
which is confirmed by the good agreement of the slope 
with Eq. ([T^ : 2) at moderate and small values of /3~^ 
the "simplified" Langevin description of Eq. ([TU]) is not 
expected to hold, and indeed discrepancy is found be- 
tween experimental slopes of 4'{A9) and those predicted 
by Eq. ((T^ : nevertheless 3) such experimental values of 
the slope appear to depend only weakly upon /3~^, so that 
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FIG. 6. Color online; Slopes of the asymmetry functions for 
several different experiments, as a function of At. Tiie inset 
shows the value of the slope (the values of each experiments 
are averaged over the plateau visible in the main plot at large 
At) as a function of /3~^ . The squares joined by the green 
dashed line represent the result of formula p3[) . 



they do not differ too much from the values at large l3~^. 
The last observation is an empirical fact which has not a 
simple explanation: at small values of /3~^ the noise felt 
by the rotator is discontinuous and the average drift can 
hardly be described as the effect of a continuous torque 
(as it is Tmotor in the FC limit). Therefore it is not clear 
at all how to define a work or an injected power and, 
consequently, a candidate for the entropy production. A 
theory for the fiuctuations of A0 in such a situation is, 
up to our knowledge, unknown and the discovery of the 
validity of the FR with a slope similar to a very different 
regime is largely unanticipated. We mention that in the 
Gaussian approximation, i.e. assuming a parabolic form 
for the large deviation rate of A9 or equivalently 



f{A6) ^ exp 



(Ae*- (to) At) 



21 



2DAt 



(17) 



leads to the identification s = 2{uj)/D. Again, no theo- 
retical expectations exist, for the stochastic process mod- 
eled in Eq. ([3]), for the ratio between the average drift 



and the angular "diffusion" coefficient D. The empirical 
observation that such a ratio is somehow independent 
from the relative importance between collisions and dry 
friction (controlled by /3~^) is quite an interesting fact. 
Interestingly, in a previous paper |23| where similar mea- 
surements have been done on the different setup cited 
in the introduction [31, the validity of the FR for the 
asymmetry function 4'{Ad) was observed at very differ- 
ent shaking strengths, with a slope ~ 0.2 independent 
from the dynamical regime (see Fig. 2a and Fig. 3b of 
that paper) and amazingly close to the slope measured 
in our experiment. 

VI. CONCLUSIONS 

In summary we have repeated the studies recently ap- 
peared in |8| and [23| , concerning the experimental mea- 
surement of the average Brownian motor effect and the 
analysis of the FR respectively, for a new setup which ap- 
pears to be simple enough for a reasonable comparison 
against kinetic theory. Such a theory predicts two main 
regimes, rare collisions and frequent collisions, with two 
different formula for the average angular velocity of the 
rotator, formula ([Ti]) and dTTT) respectively. The same 
theory is able, only in the FC, to predict the validity 
of the FR 4>{Ad) = sA9 with an analytical formula for 
the slope s, which is in good agreement with the ex- 
periments in that regime. The interesting observation, 
detailed above, is that such a slope does not depend 
strongly upon f3~^, giving a similar value even in a regime 
where collisions are rare excitations followed by fast dis- 
sipation due to dry friction. Future investigation of this 
puzzling observation is in order, in particular through 
numerical simulations or further variations of the exper- 
imental setup. 
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